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Abstract 



In this thesis we consider 't Hooft's polygon model for 2+ ID gravity. We first 
recall the ADM formalism to write general relativity in Hamiltonian form. With 
this background in mind, we give a detailed review of the polygon model and its 
explicit evolution description in the classical context. 

Then we review some remarks in the literature about quantization of this model 
and the discreteness of space-time. We discuss the problems associated with this 
in the context of canonical quantization for some explicit quantization schemes: 
the triangle inequalities conflict with the Stone- Von Neumann uniqueness theo- 
rem when we use the canonical variables. Also, the implementation of transitions, 
the constraints and their algebra poses significant problems. We conchide that 
no rigorous conclusions about the spectrum of space-time can be drawn without 
an explicit quantization scheme. 

Furthermore, we consider the Poisson structure of the constraints, which are 
important when we try to quantize the model. We improve known results and 
show that the fiill Poisson structure can be calculated explicitly and that it closes 
on shell. An attempt is made to interpret the gauge orbits generated by the 
constraints. 
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Introduction 



In this thesis we will study a particular model for (Einstein) gravity: the 't Hooft 
polygon model for gravity in 2 + 1-dimensions. We will specifically be concerned 
with the question of whether this model can be quantized and we calculate the 
complete Poisson structure of the constraints. 

The 2 + 1 here stands for 2 spatial dimensions and 1 time dimension. This 
is one dimension lower than the 3 + 1-dimensions we (seem to) live in. That 
turns out to make the whole problem a lot simpler than trying to tackle the full 
problem of four-dimensional quantum gravity, which people are trying to solve 
for decades already. It is hoped then, that by studying and hopefully solving 
three-dimensional (and in general lower dimensional) quantum gravity, we can 
learn something about how to solve four-dimensional quantum gravity. 

1.1 Gravity: general relativity 

When we say 'gravity' we mean the classical theory of gravity as formulated by 
Albert Einstein, which is known as 'general relativity'. It has as a starting point 
that (locally) the effect of a gravitational field and acceleration are the same. 

This has as consequences that space and time together (often referred to as 

space-time) can be curved: the concept of a straight line is not entirely lost, but 
becomes somewhat more subtle. For example in a curved space it can be possible 
to always travel in a straight line and come back to where you started, or that 
two straight lines intersect each other at more than one point. 

The mathematical description of this theory is much more complicated than that 
of the classical, Newtonian theory of gravity and finding solutions is much more 
difficult, already in the classical (non-quantum) theory. 

On the other hand, general relativity describes some interesting phenomena. 
Some of these are already well confirmed by measurements, like the bending of 
light rays by large masses (e.g. the sun) and the slowing of time in gravitational 
wells. Other phenomena are predicted but not yet detected, like black holes and 
gravitational waves. 
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1.2 Quantization 



At the start of the 20th century, physical phenomena were discovered which could 
not be described by the classical physical theories. Examples are the photo- 
electric effect and the observed structure of atoms. 

A theory which could correctly describe all of these phenomena was constructed 
over the years and is generically called 'quantum mechanics'. In this theory, 
observables like particle position or energy do not have definite values anymore; 
these observables have a probability-like amplitude. The 'quantum' refers to the 
fact that certain observables can only take (a superposition of) discrete, quantized 
values in this theory. 

At this moment, a theory for three out of four known fundamental forces in nature 
has been formulated: the electro-magnetic force and the weak and strong nuclear 
forces. Only the force of gravity remains "unquantized" as of this day. Yet 
one cannot consistently combine classical and quantum theories, so the current 
theories of general relativity and quantum mechanics are not complete. Hence 
physicists are looking for a quantized version of the theory of gravity. For more 
details on the subject of quantization, we refer to chapter |^. 



1.3 The polygon model 

In 3 dimensions, for solutions of the Einstein equations without cosmological 
constant, curvature of space-time does only occur locally at points where there 
is mass. Everywhere where space is empty, it is (locally) fiat. This doesn't mean 
however that things become trivial, because we can still have global non-trivial 
curvature effects. 

The polygon model for 2 + 1-dimensional gravity makes clever use of this fact. 
We can model the spatial part of space-time as a set of polygons glued together. 
Inside each polygon space is simply fiat, but these polygons can be glued together 
to form a two-dimensional spatial surface with non-vanishing two-dimensional 
curvature and non-trivial topology. We could for example create a cube from six 
four-sided polygons. The polygon model is explained in more detail in chapter ^ 

There are other models for 2 + 1-dimensional gravity. Two formulations that are 
closely related to each other and the polygon formulation are the second order 



ADM formulation with York time by Moncrief et al. [0, and the first order 
Chern-Simons formulation by Witten et al. [0, §]. On a classical level, these 
formulations are equivalent. The quantizations of these formulations, as far as 
they exist, yield different results though. Witten's formulation gives a 'frozen 
time' picture and Moncrief's formulation can only be explicitly quantized in the 
genus 1 case. 
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Notation and conventions 



In the polygon model, the fundamental variables are the lengths and boost 
parameters of the edges between polygons. These polygons cover a spatial slice 
of space-time. Boost parameters are the generators of Lorentz transformations 
between different polygons. Explicitly a boost rj is related to the speed v of a 
Lorentz transformation by 

tanh(?7) = v, 

cosh(?7) = 7(f), (2.1) 
sinh(r7) = 7(t>) v, 

where as usual we have j{v) = — - and we normalize c = 1. 

The boost parameters of edges can be viewed as Lorentz transformations in the 
specific direction perpendicular to that edge. We can more generally consider 
boosts in all directions together with spatial rotations. These form the (Lie) 
group of Lorentz transformations: we can compose any two elements to get 
another Lorentz transformation. We write 0(2, 1) for the group of Lorentz trans- 
formations in three-dimensional space-time. This group can be generated^ by the 
three basic transformations of a rotation and boosts in the x and y-direction: 



LJa) 



I coshr] sinh?/ . 

Lxiv) = I sinh?] coshr/ , (2.2) 
1 

cosh 1] sinh r] 

Lyiv) = I 1 

sinh rj cosh 77 




^Actually this will only generate the identity component 5*0+ (2, 1) C 0{2, 1), but this will 
be sufficient for our purposes; see appendix ^ for more details. 
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As generators of Lorentz transformations, boosts have the nice property of addi- 
tivity: 

Lxivi) Lx{V2) = L^ivi + V2)- (2.3) 

We will denote boosts by Greek letters from the middle of the alphabet {i], fi,u, . . .), 
angles by Greek letters from the start (a, 7, . . . ) and lengths by l,m,n, . . . 

In the polygon model, we will frequently encounter hyperbolic functions of the 
boosts and trigonometric functions of angles. To keep formulas short and clear, 
the following notation will be used: 



Vi 




= tanh(?7j). 


li 


= Im 


= cosh (2 rji 




= 


= sinh(2 


Ci 


— c-di 


= cos(aj). 


Si 




= sin(ai). 



Note the extra factor 2 in the definition of 7 and a. Furthermore, when the 
specification of a variable is obvious from the subscript index, then this shorter 
notation will be used. 

Boosts, angles and lengths will sometimes be indexed along the border of a 
polygon and sometimes along a vertex. To distinguish between those, we will 
label the former with indices i,j,k and the latter with indices a,b,c. Boosts 
and angles with the same index are related in different ways for those cases. 
At a vertex, a boost and its opposite angle have the same index, while along a 
polygon an angle and the boost of the counterclockwise following edge have the 
same index. See figures |2.1j and |2.2| . Note that when edges and/or vertices have 
to be numbered, this will always be in a counterclockwise orientation. 




In, Vn 

Figure 2.1: Labels of edges Figure 2.2: Labels of edges and angles 

and angles at a vertex. along a polygon. 
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Furthermore we see that indices along a polygon are cyclic. In computations we 
will therefore always implicitly interpret the index numbers modulo (A^ being 
the number of vertices and edges of the polygon). For example, we have that 
Vo = Vn- An exception to this rule will be made for quantities that are sums 
over ranges of indices, as a sum over the empty range is truly different from a 
sum over the full cycle. For example, the sum of angles up to index j as defined 
(4.6)| : here we interpret 



m 



^TT-a, = (2.5) 
1=1 



instead of 6*0 = 6* 



N- 



The model has a phase space consisting of a set of edge lengths and boosts (/j, rjj). 
The Poisson brackets are defined on this phase space as 

{k,r^^} = l6,,, (2.6) 

where the indices i,j run over all edges. 

The phase space of this model is larger than the physical phase space. This is 
reflected by the presence of algebraic constraints on the phase space variables. 
These define a subspace of the complete phase space by imposing C{l,r]) =0 for 
all constraints C. When dealing with constraints, we use the '~' sign to indicate 
an equality, that holds only on this constraint surface and not necessarily on the 
whole phase space. Thus by definition a constraint function 0(1,1]) satisfies 

C{l,v)^0. (2.7) 

The polygons that together make up a 'tessellation', can be viewed as a graph 
structure, consisting of faces (the polygons), edges (the polygon edges) and 
vertices (the points where three or more edges join). The sets of these faces, 
edges and vertices will be denoted by F, E and V respectively. Single elements 
we denote by / G F and the number of elements by See section [A.l| for 

some more details. 



A note on signs 

Some of the definitions used in this thesis differ by a minus sign from referenced 
works or conventions; unfortunately it is not possible to choose definitions to be 
simultaneously convenient, according to conventions and equal to all referenced 
works. To keep things clear, we give a (hopefully complete) list of differences. 
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The Poisson bracket has a relative minus sign with respect to the definition in 
[Q; it is chosen to have the same sign as the standard definition. Related to that 
is the definition of the time evolution of a classical quantity (function on phase 
space) as 

f = {f,H}, (2.8) 

which is also according to normal conventions and opposite to the definition used 
in p and related papers. 

Boosts have the same sign as in contracting edges have positive boosts. 
This looks a bit like a counter-intuitive choice, but is necessary to get the right 
Hamiltonian fiow for our choice of signs in (2.6) , K2.8) and the Hamiltonian. 

We use the space-like (— , +, +) sign convention for the metric. 
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Representing gravity 



Before we can start to think about the quantization of gravity in arbitrary 
dimensions, we first have to choose a way to represent gravity. We would hke to 
be able to formulate gravity within the framework of Lagrangian or Hamiltonian 
mechanics. Or in other words, we want choose a set of variables that fully describe 
a state of the system and then formulate the equations of motion on those states. 

We must note however, that a classical "state" is normally defined as a full 
description of the system for a given instance of time. But in the context of general 
relativity "a given instance of time" is a priori not an unambiguous statement to 
make: there is no unique time coordinate, because coordinate diffeomorphisms 
are a symmetry of the system. 

In the next section we will introduce a way to split time and space in general 
relativity, which then makes it possible to formulate general relativity in a La- 
grangian and then Hamiltonian formulation of gravity. 

3.1 Foliation of space-time 

Before trying to split off time, we first have to give a mathematical description 
of a space-time. There are two well known choices of fundamental variables to 
describe a general space-time. We will use the standard choice of taking the 
metric g^y{x^) as field variable. There is however an alternative formulation, 
known as 'first order formalism' in which a set of orthonormal vectors (called 
'dreibein', 'vierbein' or in arbitrary dimensions 'vielbein') together with an in- 
ternal connection is chosen as fundamental variables. This description relates to 
the standard metric description via 

r]ahelel = g^^. (3.1) 

In space-time dimension n, rjab is a diagonal matrix with signature (1, n— 1), which 
expresses the fact that the are a local orthonormal set of vectors. 

A space-time can mathematically be specified in a very general manner as a 
differentiable manifold M with a metric g on it, often denoted as a (pseudo-) 
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Riemannian manifold (M, g). In the following, we will be dealing with a pseudo- 
Riemannian manifold of signature (1,2). 

We now want to formulate the theory in a Hamil- 
tonian formalism, which can later be used 

starting point for canonical quantization. For this, V'^^'^^^^-Cj' ~ -1 - _ 

we will choose an explicit time coordinate and " ^^fTT — -3 

introduce the concept of a global foliation of space- " ~ - - f 

time. This can be seen as a continuous set of slices ~ ~ -^^TT — [ 

of space-time, where the chosen time parameter 

labels the different slices; time has a fixed value Figure 3.1: A foliation: the 

on the whole of each slice, see figure O. To this indicates the direction of 



time. 



end, we introduce a time function T on M, asso( 
ating with each point in M a unique time t. This 

function should be sufficiently well-behaved such that its original = T~^{t) for 
each instant of time is a submanifold of M, which is spatial. It must also be a 
monotonically increasing function along all timelike curves. 

The process above amounts to a partial choice of space-time^ coordinates: we 
have chosen (only) a time coordinate, which slices the manifold into a set of 
submanifolds of equal time. Note however that this "time" is not necessarily the 
proper time for a local observer! This complete set of submanifolds we now call 
a 'time foliation' of the manifold M. Note that there is no unique foliation: any 
time function T that satisfies some regularity properties will suffice and give a 
different foliation. 

Now that we have a time foliation, we want to pick a set of configuration variables 
that for a given time t specify the state of the system, just like a function x{t) 
would describe the state (position) of a pointlike particle for a given time t. 
Heuristically we can say that the system is the metric on the space-time manifold 
M and at a given instant of time t, it is given by a spatial slice Sj with the 
restricted metric. The configuration variables for the complete manifold are given 
by gf_iu and for a fixed time t by restricting the metric to and only looking at 
the space-space part of the metric: that is the part needed together with to 
describe the state of space-time at time t. We also need to specify the conjugate 
momenta to get a complete set of initial conditions. 

When we choose a complete set of space-time coordinates, we can make this 
choice of configuration variables explicit by defining 

9,u=\ f \\- (3-2) 

^One must be careful to make the distinction between coordinates that map space-time and 
coordinates that describe the state of the gravity system (which we will refer to as 'configuration 
variables'). 
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The A^, N"^ and Qij are functions of space-time that can be read off from the 
exphcit form of g^^, in the chosen space-time coordinates. We have chosen time 
to be a global timelike coordinate function, so this implies that goo < and 
also that the symmetric submatrix gij for the other spatial coordinates is non- 
degenerate and the entries gij thus well-defined. 

The function is in the literature (see e.g. [@, p. 12]) referred to as the 
'lapse function': it specifies the rate that clocks tick with respect to the chosen 
time coordinate. The functions A^* are called the 'shift functions' and specify 
the displacement of spatial coordinates, when we make a small displacement 
perpendicular to the surface Sf. The matrix gij can be viewed as a metric on the 
submanifold and thus describes the state of the system (to be compared to 
the position x in a one-particle classical mechanical system). 

Thus far, we have not yet made any specific choice of space-time coordinates, 
except for the fact that we demanded a global time coordinate. The only thing we 
have done is to make an explicit separation of the time from the space coordinates 
to allow for a treatment of the theory within a Lagrangian and Hamiltonian 
framework. This decomposition is known as the Arnowitt-Deser-Misner (ADM) 
formalism. 



3.2 Hamiltonian formulation 

Using the ADM decomposition, we can write down the action of general relativity 
and transform it to a Hamiltonian formulation. We will here give an outline of 
this procedure, to give some insight in the intricacies that will also show up in 
the polygon model. More details can be found e.g. in reference [0]. 

The action of general relativity in three dimensions is given by 

^ = 16^ / "^'"^ " ^^-^^ 

in the absence of external matter or fields. We normalize IQnG = 1 and set the 
cosmological constant A = from here on. 

Inserting the ADM decomposition of the metric, this action can be rewritten as 
an integral over time t of a Lagrangian L, with 



/ d^xN^\\R + K,^K'^ - [g'^K.^Y] . (3.4) 



Here R is the Ricci scalar corresponding to g and Kij is the extrinsic curvature of 
Ht embedded in M. In principle there are also some boundary terms from partial 
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integration, but these cancel because the spatial slice is compact. The integrand 
3.4)1 we denote by £, the Lagrangian density. 



m 



Now we have a standard Lagrangian system and can make a transformation to 
canonical momenta, which yields 



dC 



with 



/ 

n 



d [dtgij] 



dt / d\ {tt'^ dtgij - NH- NiH') 

St 



-2 Vj ir'^ 



~g\R, 



The complete Hamiltonian, given by 



H = [ d^x {Nn + NiW) , 



(3.5) 

(3.6) 

(3.7) 
(3.8) 

(3.9) 



is now a linear combination of constraints, because the functions A^, A^* do not 
have any equations of motion associated with them and thus act as Lagrange 



multipliers for 7i,7i*: the variation of the action |(3.7)| with respect to N,N^ 
yields H = 0, W = 0. 



3.3 Reduction of phase space 

We have found that in the Hamiltonian formulation we are left with constraint 
functions. This means that the real dynamics of the system takes place only in 
a subspace of the full phase space that we had chosen initially. 

We started with a configuration space Q of all metrics cjij on a spatial slice and 
the corresponding phase space of canonically conjugate pairs (^jj,7r^'). To get a 
description of only the true physical degrees of freedom, we have to reduce this 
phase space. This is done by first solving the constraint functions |(3.7)| and |(3.8) 



In addition, the constraints above are first class constraints, so they generate 
gauge transformations. The functions Ti^ generate coordinate transformations of 
the surface S^, which is a well known gauge freedom of general relativity. The 
function Ti generates translations in time, relating coordinates on different time 
slices S. This is in the strict sense not a true gauge transformation, because 
it relates configuration variables at different times. This fact is related to the 
problem of time in (the quantization of) general relativity [|l^ . 
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Gauge transformations do not change the physical state of the system, so as a 
second step, we have to fix this gauge freedom to get the reduced phase space. 
The gauge freedom generated by the first class constraints (3.7) and |(3.8)| can now 



be removed by fixing the Lagrange multiplier functions that were still arbitrary. 

The explicit work of solving the constraints and fixing the gauge is not as straight- 
forward as it might seem from the formulation above: there is no known way to 
solve this in 3+1-dimensions, but in 2+1-dimensions there is [p!^ , p!9[| . We will not 
go into the details, but with the use of Riemann surface theory one can consider 
the space of metrics g modulo diffeomorphisms and conformal factors, which is 
a finite dimensional space known as the moduli space A/" of S. This space has 
dimension 6{g—l) when S has genus g > 1, two when g = 1 and zero when g = 
(see also figure for the concept of genus of a surface). The physical phase 
space now has twice the dimension of this moduli space A/'(S), because we still 
have to add the canonical momenta, which span the same number of dimensions. 

Thus in 2 + 1-dimensions, the physical phase space of general relativity is finite 
dimensional, unlike the phase space we started with: that was the space of all 
metrics g and conjugate momenta, which are fields over E and thus infinite 
dimensional. For this reason, general relativity in 2 + 1-dimensions is sometimes 
called a 'topological theory': the true degrees of freedom do only arise from 
topological non-triviality of the surface S. This again is a reason that the polygon 
model works: there are no local degrees of freedom and locally space-time is fiat. 
The formulation of the polygon model depends on this fact as we will see in 
chapter |^. 
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The polygon model 



The polygon model is a model to describe 2+1-dimensional space-times that can 
be foliated by a set of spacelike Cauchy surfaces. This includes the possibility of 
spatially open and closed universes although in the case of an open universe 
one has to be careful about boundary conditions at spatial infinity. 

We will however only consider spatially closed universes. Thus if we take a spatial 
slice of our complete space-time manifold M, then this is a closed and bounded 
surface (so it has finite volume). More specifically, we take a space-time of the 
form of a product: 

M = S X /. (4.1) 

Here / C R is some interval of time, which may be bounded from below and/or 
above, depending on whether there is a big bang or big crunch respectively. 

S is the spatial part and thus should be a compact two-dimensional surface. These 
surfaces are topologically completely classified. The ones that are orientable are 
classified by their genus g, which indicates the number of 'holes' in them. The 
simplest examples are a sphere, a torus and a 'double torus', see figure 

The model also allows for a finite number of pointlike, spinless particles to be 
added to this space-time. We have not studied this in detail and refer the reader 
to p, |lT[ for more information. 

The simplest spatial surface of a sphere is only realizable with particles. Proof 
for this fact can be found in p. 57] or [0, p. 2912] and is based on an analysis 
of the moduli space of metrics on the surface S modulo diffeomorphisms. For the 
sphere this moduli space consists of one unique metric with curvature 1, which 
does not yield a solution to the vacuum Einstein equations. This can also be seen 




Figure 4.1: Surfaces of genus 0, 1 and 2. 
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from the fact that the Hamiltonian constraint can be written as in |(4.ir 

^ = 2^X = 47r(l-(7). 



The Hamiltonian equals minus the area of the dual graph as given in |(A.10' 



This area cannot be negative, thus the Hamiltonian cannot be positive, which 
rules out the case g = 0. 



4.1 Tessellation of the spatial surface 

We are considering general relativity in 3 dimensions. This we can exploit to 
simplify the description of our space-time. In this case the Riemann tensor is 
completely determined by the Ricci tensor: one can show this e.g. by counting 
the degrees of freedom or explicitly p, p. 3]: 

Rpupa 9py Rua ~\~ 9ua Rpp 9yp Rpu 9pu Rup 2 ^9pp 9^^ 9pcr 9vp) R- l^'^) 

This means that in a region of space-time with no matter (T^'^ = 0), the Einstein 
equations effectively reduce to R^^ = and determine that space-time is locally 
flat. Thus everywhere we can locally choose a Minkowski coordinate system. 

Given this space which is locally Minkowski, we can also easily choose a spatial 
slice that is locally fiat too. This will induce a local coordinate system from the 
three-dimensional Minkowski space to our two-dimensional spatial slice S^. We 
cannot in general extend this coordinate system to our whole spatial surface: 
our space is compact, so when we try to extend our flat coordinate system from 
a certain starting point, we will run into parts of space where we had already 
defined coordinates. These do not necessarily have to match; while we have 
extended our coordinates and run into a "charted" part of space, we might have 
gone around a non-contractible loop in S which has a non-trivial holonomy. A 
simple example is a toroidal spatial surface which expands uniformly in one of its 
fundamental directions: when we walk along a loop in this direction and return 
to our starting point, the coordinates fail to match by a Lorentz boost of the 
speed of contraction or expansion. 

We can however choose a coordinate system at some point (s) of our spatial 
surface and extend it as far as possible. After doing this, we will be left with 
several regions of fiat space which, glued together at their boundaries, form 
the complete spatial surface again. We will call such a set of fiat coordinate 
patches a 'tessellation'. Notice that at points where regions of space are glued 
together, the spatial surface will in general not be fiat. The underlying three- 
dimensional Minkowski space-time is fiat, but at points where three or more 
spacelike slices join, the two-dimensional surface will in general have a curvature 
of delta distribution type. See figure for an example. 
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Figure 4.2: Curvature at a vertex 
of three joined 2D surfaces in a flat 
3D background. 



Figure 4.3: Tessellation of a torus 
by 2 polygons, labelled I and II. 



This is a tessellation of a spatial surface at one fixed value of the time coordinate. 
The coordinate systems of the different patches should therefore obey some 
boundary conditions, such that there is no time difference when traversing a 
boundary and also time should run equally fast on both sides. This implies 
that the boundary should move perpendicular to its orientation in the coordinate 
frames on both sides, which yields that a boundary must be a straight line and 
its velocity be of same size and opposite direction in the neighboring coordinate 
frames (see [|, p. 1336] for a more detailed discussion). 

An example of a tessellation is given in figure the dashed lines show a 
fundamental domain of the torus (with the usual identification of opposite sides) 
and two 6-sided polygons covering this torus. The two polygons are labelled I 
and II and the different edges are labelled 1 to 4. Notice that these polygons 
have a number of edges in common with each other and also with themselves. 

Thus our tessellation consists of a number of polygon coordinate patches with 
sides that have to be identified. The sides of these polygons have well-defined 
lengths within the coordinate frames of the respective polygons. These should 
however match the length of the edge as seen in the coordinate frame on the op- 
posite side of the boundary (because the metric should be continuous). Therefore 
these boundary edges have a unique, well-defined length /. 

Besides that, we can assign a boost parameter to each edge. This boost is = 
tanh~^(t;), where v is the speed of the edge in its neighboring coordinate frames. 
The speed and boost have a positive sign when the edge moves inward, thus v 
is the speed at which the edge contracts the polygon. This speed is uniquely 
defined, because an edge must expand or contract with the same speed in the 
coordinate frames on both its sides. The Lorentz transformation relating two 
neighboring coordinate frames is then given by a boost of 2 rj: first a boost from 
one frame to the rest frame of the edge and then again a boost from the edge 
rest frame to the other coordinate frame. 
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4.2 Vertex relations 



We have made a tessellation of our spatial slice with polygons. These polygons 
exactly cover this spatial slice, so boundaries ('edges') of these polygons have to 
be identified where they are "glued together" . 

Besides that we also have 'vertices': points where three or more edges join. As 
was stated earlier, space-time is locally fiat. This should also hold at vertices. To 
express the curvature located at a vertex, we can write down the holonomy around 
that vertex: the change a vector undergoes when it is being parallel transported 
around that vertex. 

We can explicitly write down the complete holonomy of going around the vertex 
using the Lorentz matrices |(2.2)| . Going around in anti-clockwise orientation. 



we label the polygons at the vertex i = 1 . . .n and define aj as the angle that 
polygon i makes at the vertex and 2'r]i the boost between polygon i and polygon 
i + 1 mod n. Now going around means successively making a rotation Lr{ai) and 
a boost Lx{2r]i) for each polygon.^ This should be equal to the identity operation, 
as space-time must be fiat at the vertex point, thus we have 



i=l 



'a,)L^{2rj,) = 1. (4.3) 



We demand that at a vertex always exactly 3 polygons join. This is always 
possible by choosing a suitable tessellation and will allow us to express the angles 
in terms of the boosts. For later convenience, we now choose the indices of boosts 
and angles in a more symmetric fashion as in figure |2.1| . When we rewrite |(4.3) 



to a form with 3 matrices on the left and the right side and then explicitly work 
out the equations, we get: 

si : S2 : S3 = CTi : 0-2 : as, (4.4a) 

= S1C2 + 72S3 + C1S273, (4.4b) 

ci = C2C3 - 71S2S3, (4.4c) 

7i = 72 73 + ci a2 0-3 (4.4d) 

and all cyclic permutations of the indices. Here we used the notation as in 



2.4)[ As in the literature, we will refer to these as the 'vertex relations'. These 



equations are not independent nor have we given the same set as originally 

^Remember that the transformation between two neighboring polygons goes with twice the 
boost of the edge. Furthermore the precise choice of signs and Lx or Ly depends on the exphcit 
choice of coordinates in the polygons and orientation of the holonomy, but this does not affect 
the resulting equations. This can be seen from the fact that the system has space and time 
mirror symmetry. 
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given by 't Hooft in p , but they are complete in the sense that they are equivalent 
to the full set of equations determined by |(4.3) 



These equations now give us expressions for the angles of polygons in terms of 
boost parameters or vice versa. To determine the angles one can for example 
rewrite equation |(4.4d)| and take the inverse cosine. This does not yet determine 
the angle uniquely, but with the help of equation |(4.4a)| it does. There is the 
extra fact that at most one of the three angles at a vertex can be larger than vr, 
because if there would be two angles larger than tt, there would be points that 
show up in the coordinate systems of both polygons This then determines 
the sign of the sines of the angles at each vertex. 



4.3 Constraints 

We now have a representation of a 2D spatial slice of the 2+ ID space-time 
manifold M in terms of one or more polygons, glued together at their edges. 
A configuration of such a spatial slice is described in terms of the lengths k and 
boost parameters rji of all edges together with the graph structure of all these 
edges. 

The li,rii however cannot take all possible values in M^^. Besides the fact that 
all lengths have to be positive {li > 0), there are some other, more complicated 
constraints. 

Firstly the boosts at each vertex have to obey a triangle inequality. For each 
triple of boosts at one vertex and each permutation thereof, the following must 
hold: 

\Va\ + \Vb\ > \Vc\ (4.5) 
This can be deduced from vertex relation |(4.4d)| . A more geometrical argument 



is as follows: when going around the vertex, we are Lorentz boosted three times, 
but we should return to our original rest frame. The resulting boost of two boosts 
is at most the sum of the boosts, exactly when they are coUinear. So the sum 
of lengths of each two boosts must be greater than the third to be able to close 
the holonomy. This is analogous to the case of Galilean transformations when we 
replace the boosts by normal velocities. 

Secondly there is a set of constraints for each polygon. Geometrically these are 
easy to formulate: each polygon must close. This means that the angles and 
edge lengths must be such that when going around the border of a polygon and 
keeping track of the relative coordinates, that one must end up exactly where 
one started and the sum of the angles must be 2 n. 

As the angles can be expressed in terms of the boosts, we can explicitly express 
these constraints in terms of the model variables li,rii. For a fixed polygon, let 
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us number the lengths, boosts and angles according to figure 

We first define the angle between edge j and the horizontal (edge A^) by 



9, = J2^-a.. (4.6) 



1=1 

Using this we can write down the angular and closure constraints as 

N 



Ce = On -2n = ^(vr -aj)-2n^ 0, (4.7a) 

N 

C, = J2 h exp(« Oj) ^ 0. (4.7b) 



The closure constraint is written using an implicit complex coordinate system 
in the polygon. This way of writing is obviously equivalent to formulating two 
real-valued constraints for the x and y direction. We will refer to the second 
constraint as the 'complex constraint' as in [ffH]. 



4.4 Dynamics 

From the definition of our length and boost parameters, we can find the evolution 
of these parameters. The boosts do not change in time, except when transitions 
take place. The lengths however do change. A geometrical analysis shows that 
the length change of an edge gets a contribution from the vertices at both its 
sides equal to 

dl tanh(?7i) cos(a3) + tanh(772) ''^i cs + ^2 / , 

at sm(a3j S3 

where index 1 numbers the edge under consideration and indices 2 and 3 may be 
interchanged (there is mirror symmetry). 

We can also look at the evolution of the system in the Hamiltonian formulation. 
The Hamiltonian is given by the total two-dimensional curvature of the spatial 
slice. Curvature is only present at vertices and is given there by the deficit angle, 
thus 

^^ = 5Z (27r-^a,,| . (4.9) 

We can see that this Hamiltonian is a constraint: on the one hand, it is the integral 
over the scalar curvature, which is only dependent on the Euler characteristic x 
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(see formula 



(A. 2)1 ) by the Gauss-Bonnet theorem: 



H 



j^d'x ^R = 27rx = ^7r{l-g). (4.10) 



On the other hand we can see that it is equal to the sum of all angular constraints, 
using the Euler characteristic of the graph: 



= E(^/-2 
= -^a, + (3#V-2#F)7r 



= E 2vr-5^a., +(#y-2#F)7r 

= H-27rx. (4.11) 

In the last step we made use of the trivalence of the vertices by plugging in the 
relation 3 = 2 ^E. See page ^ for notation. 

The Hamiltonian now generates evolution of the system by 

§ = {k,H}, 

(4.12) 

-^ = {V.,H}. 

To explicitly calculate this, we need a symplectic structure on phase space. A 
symplectic structure defines Poisson brackets and it is uniquely determined by 
the Poisson brackets between the fundamental variables. In our case it turns 
out that the evolution of the boosts and lengths exactly matches the evolution 



generated by the Hamiltonian if we choose the symplectic structure to be |(2.6 



The Hamiltonian time evolution of the boost parameters can immediately be 
seen to match that of the geometrical picture: in both cases it is zero. For 
the geometrical picture it was chosen this way by letting time run equally fast 
everywhere in each polygon. For the Hamiltonian evolution, we see that it follows 
from the fact that H only depends on the boost parameters. 
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4.5 Transitions 



Given the results of the previous section, the evolution of the system looks simple: 
the boosts and angles are fixed and the lengths change linearly in time. It turns 
out that things are not that simple: there is a complicating factor coming from 
the constraint that the model parameters must describe a true tessellation. This 
means that edge lengths must be positive and polygons must be true, non-self- 
intersecting polygons. 

As a side remark, it must be noted that strictly speaking, polygons may be 
self-intersecting, but only in such a way that they can be embedded in a non-self- 
intersecting way on some two-dimensional surface, see P]. This can be formulated 
as follows: if a polygon intersects itself then this will generate a transition if and 
only if a new triangle is formed and the winding number of points that will lie 
in the triangle decreases with respect to the oriented boundary of the polygon. 
If the winding number of these points decreases, then these points are on the 
outside of both parts of the polygon boundary that cross and thus would start 
to show up in two different places in polygon coordinate charts. 



Figures [4.4| and |4.5| show two polygons that both seem to self- intersect. In the 
first a transition will split the polygon into two parts when the vertex hits the 
horizontal edge. The winding number of points in the new triangle decreases 
from to —1. In the second figure, no transition will take place when the vertex 
crosses the edges as indicated. 

Returning from this digression, we now have essentially two transitions that can 
happen: an edge collapsing to zero length (the top left diagram in figure and 
a vertex hitting an edge (top right). There are some special cases of these two 
transitions. If an edge is part of a triangular polygon, then this whole polygon 
will disappear when the edge length goes to zero. This is due to the fact that 
boosts and angles are constant and thus a triangle will scale to zero size when 
one of its edges does. This can also happen with two neighboring triangles (but 





Figure 4.4: A truly self- 
intersecting polygon. 




Figure 4.5: A superficially 
self-intersecting polygon. 
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Figure 4.6: All possible polygon transitions without particles, labelled as in j^. 



not with three, as they would together close the spatial slice, but violate initial 
conditions). See the middle right and bottom diagrams of figure [4 .61. Furthermore 



we have a special case when an edge shrinks to zero, but the angles are such that 
it starts growing again, but with angles changed by tt (called a 'grazing transition' 
in the literature); see the middle left diagram of figure ^]6|. 

There are a few more possible transitions when one allows particles in the model. 
For the full set of transitions including those with particles we refer to 0. 

Now if we have a set of initial conditions that satisfy the constraints, we can 
calculate the time evolution. This is in first instance given by only the length 
change of the edges. But after a certain amount of time, the polygon configuration 



may have changed in such a way that a transition as depicted in figure ^l6| is about 
to take place. This will not only change the graph structure as indicated, but 
also the model parameters will change. One can find these by demanding that 
the parameters after the transition should be consistent with those before and 
that they should obey the vertex relations: space-time must stay fiat. 

A potential problem now arises: when a transition takes place, we switch to a 
different graph structure with a different set of edges and thus a different phase 
space to describe the system. The number of edges does not need to be the 
same, so the dimension of the phase space can change during a transition. The 
dimension of the physical phase space does not change (as expected). This we can 
see from counting degrees of freedom and constraints. The Euler characteristic 



(see section [A.l|) is 

X(S) = #F - #i? + #F (4.13) 
and is constant for a given topology of S. From the trivalence of the vertices, we 



22 



obtain 3 = 2 ^E. This together gives us the relation 

#F - 3 #^ = constant, (4.14) 

which does not depend on the specific structure of the graph. Now two degrees 
of freedom are associated with each edge. Each polygon on the other hand 
introduces three constraints. And as we have seen in section p.3| , each constraint 
generates a complementary gauge degree of freedom. So each polygon amounts 
to a total of six unphysical degrees of freedom in the full phase space, which 
exactly cancels the extra degrees of freedom added by the edges. 



4.6 Constrained dynamics 

With the addition of transitions, we now have a complete description of the 
(classical) dynamics of the system: given a set of initial conditions, we can 
uniquely determine the evolution of the system. This system must not only 
be uniquely defined, but also be well-defined. For this we have to check that all 
constraints (equalities and inequalities) are preserved under the evolution of the 
system. 



4.6.1 Equality constraints 

We want to check that time evolution preserves the constraints, which amounts 
to checking that { C , } ~ for all constraints. This follows trivially for the 
angular constraints, because they do not depend on the lengths. The calculation 
for the complex constraints is not straightforward, but also gives { , if } ~ : 
see chapter ^ for the details. 

Besides this, one must also check that the constraints are preserved during 
transitions. This turns out to be relatively simple for the complex constraint: 
during transitions of the type where an edge length reduces to zero (transitions 
a, f, g, j in figure [4.6|) , only edges with length zero are removed and created. Thus 
these do not affect the complex constraints. Also in the case of transition type 
d, the complex constraint is conserved: the new lengths of the split edges are by 
construction chosen to match the complex constraint. 

The angular constraints are also preserved under transitions. The best way to 
understand this is by looking at the dual graph (see appendix ^ and p!3|). 
This graph is constructed by mapping vertices into faces and vice versa. Edges 
are mapped to themselves and to each edge in the dual graph we associate an 
(oriented) length 1 = 21], where r] was the boost of the edge in the original graph. 
To each angle a we associate an angle a = n — a in the dual graph as shown in 
figure ffTT. 
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Figure 4.7: An edge transition and the dual graph (dotted hnes). 



Now a transition will also induce a corresponding transition in the dual graph. 
As explained in appendix ^ the dual graph (with the mapping of edge lengths 
and angles) can be embedded in the hyperbolic plane. From the transitions in 
the dual graph, it can now be shown that the angular constraints are preserved 

Making an edge transition of type a (figure |4.6| ) is now a matter of removing 
an edge in the dual graph and replacing it by the opposite edge within the 



quadrilateral just created, as in figure 4.7. This dual graph is embedded in 



hyperbolic space, where every point has a full 2 tt neighborhood (as in fiat space). 
Each neighborhood of a vertex in the dual graph is also fully covered before 
the transition, as can be deduced from the angular constraint of the polygon 
corresponding to that vertex. 

The other transitions can be dealt with in a similar fashion: the angular con- 
straints are encoded in the dual graph as vertices containing a 2 vr neighborhood 
and making the corresponding transition in the dual graph does not change this, 
thus the transition should preserve the angular constraint. We have not checked 
this explicitly and there might arise some problems, when mixed vertices are 
involved, see section . 



4.6.2 Inequality constraints 

In section [4.3| we saw that there are also some inequalities to be satisfied, namely. 



lengths must be positive and the boost parameters must obey triangle inequal- 



ity 1(4. 5)[ This raises some problems when trying to find a quantized version of 



this model. 

In a standard Hamiltonian formulation, we have a configuration space Q and a 
state of the system is given as a point in the cotangent space T*Q. Normally Q 
is chosen to consist of the positions (lengths in our case) and then elements in 
the vector space T^Q are the momenta (the boosts in our case). We now face the 
problem that the Hamiltonian is only defined for those boosts that satisfy the 
triangle inequality, and thus H is not a function on the whole phase space T*Q. 
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Figure 4.8: The restriction of 
boosts by the triangle inequahties. 



Figure 4.9: Parameters at a 
vertex with one zero boost. 



In figure ^78] , a picture is given of the triples of boosts (?7a, ?7c) at one vertex 
that satisfy these inequahties. These are the triples within the triangles like the 
one shaded in the figure, stretching out radially. Note that the picture does only 
display the triples in the first octant of (all boosts positive), so the full picture 
should be mirrored in three planes. Moreover, it only shows the constraints from 
one vertex; each edge is connected to two vertices, making the complete set of 
inequalities more complicated. 

Special care has to be taken at the boundary of the inequalities, when they are 
exactly satisfied. Here we can distinguish two cases according to whether one of 
the three boosts rja, r]b, rjc is zero or not. If none of the three is zero, the angle 
opposite to the largest boost is zero, which does not correspond to an acceptable 
tessellation. 

If one of the boosts is zero, the other two must be equal (also in sign) and a degree 
of freedom in the boosts is replaced by a degree of freedom in the angles at that 
vertex. This can easily be seen from the geometrical picture (see figure |4.9|) : if one 
edge boost is zero, the polygons separated by that edge can be joined together as 
the Lorentz transformation between those polygons is trivial. Then the two other 
edges with boosts r^i , 772 become one edge and thus should have the same boosts 
rji = r}2 and angle vr between them. In this case, the angle a is not determined 
by the vertex relations and can be chosen freely. Given a, we have 13 = n — a. 
If we also take 771 = 0, the triangle inequality implies that also 772 = and the 
three angles are free but should add up to 2 tt. 

This change from boosts to angles to parametrize the degrees of freedom in phase 
space seems to indicate that the /i,?]^ coordinates cannot cover the full phase 
space. The subspace not covered seems to be lower-dimensional though. 
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4.7 Simulation 



Lastly, it must be mentioned that this polygon model has as one of its virtues 
that it is well suited for numerical simulation 0. The model has an explicit 
description of the evolution of the system, which is finite dimensional and not in 
terms of differential or integral equations. This allows one to "exactly" simulate 
it on a computer: the exactness is only bound by the finite precision used in the 
simulation. 

To simulate a space-time, one must first find a set of parameters that satisfy 
the constraints. This is in general very hard, but an explicit (but intricate) 
method for genus g > 1 has been given by Kadar and Loll in |jl3[ using the 
dual-graph representation (see appendix A more simple approach is to start 
with a completely symmetric one polygon tessellation (OPT) and perturb the 
parameters slightly. One can then pick one boost and two lengths and solve the 
angular and complex constraints for those. This approach does not yield initial 
conditions for the full phase space however. 

Now, one can simulate the polygon model by calculating the length changes 



li from 1(4.8)1 and let time run. If this length change will result in one of the 



transitions in figure [4.6|, one has to "stop" the simulation and perform the 



transition, by changing the graph structure and calculating the new parameters 
using the vertex relations (4.4) . This can then be repeated until we either run 



into a big crunch, or there are no transitions taking place anymore and space 
expands to infinity. One can also reverse time by changing sign of all boosts, to 
simulate the past of a given configuration. 

We mentioned that this simulation is exact up to numerical precision. Unfortu- 
nately, if the universe collapses towards a big crunch, the boosts tend to grow 
exponentially under transitions and precision is lost after only a couple of tran- 
sitions. This inaccuracy cannot be prevented easily as the boosts are calculated 
using hyperbolic sines and cosines, thus making intermediate expressions roughly 
exponentials of the boosts and absolute errors blow up equally fast. 

We have implemented a simulation program too, to gain a better understanding of 
the details of the model and with the hope of revealing some interesting aspects 
of 2 + 1-dimensional gravity for genus g > 1. This last goal has showed to be 
fairly difficult: the interpretation and visualization of simulation data is hard 
and especially the extraction of physical information turned out to be difficult. 
The observations that we made in this simulation as mentioned above, do coincide 
with those already made by 't Hooft in M. 
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Quantization 



The polygon model for 2 + 1- dimensional gravity is a relatively simple model, 
with only a finite number of degrees of freedom. Furthermore we have been able 
to explicitly formulate it in a Hamiltonian formalism. This raises the question 
whether it is possible to quantize this model. Even though a full four-dimensional 
quantum gravity theory would still be far away, it might give interesting indica- 
tions about qualitative aspects of quantization of general relativity. 

One such thing is the question whether space and/or time will become discretized. 
Already for this model there are different predictions: 't Hooft argues in that 
the imaginary part of (analytically continued) lengths will be quantized and that 
time will be discretized. This means that the time evolution is well-defined for 
discretized time steps only. Waelbroeck on the other hand argues in that this 



time discretization can be lifted by choosing an internal time: the Hamiltonian 
constraint is solved for one of the boosts and this boost, expressed in terms of the 
other boosts, will then take the role of the Hamiltonian and its conjugate edge 
length will take the role of time. 

As already argued in the papers cited above, the quantized version of this theory 
might very well depend on the way we quantize. We have a set of constraints, 
which we can choose to first solve and only then quantize the reduced phase 
space. On the other hand, we can also first quantize and afterwards impose these 
constraints as operators in the quantum theory. 

These predictions about a quantum theory of the polygon model are heavily 
dependent on the way one quantizes the theory. We will be concerned with find- 
ing a canonical quantization only and not consider other quantization methods 
like path integral quantization: canonical quantization is already a non-unique 
quantization scheme, that is also not guaranteed to work. Therefore we will first 
investigate the a priori question of whether a consistent canonical quantization 
exists, before further investigating the predictions of a quantized model. 
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5.1 General procedure 



Here, we will describe what it means in a general setting to canonically quantize 
a classical system. 

We start with a classical system with a configuration space Q, which is locally 
isomorphic to M", or in other words, which is an n-dimensional smooth manifold. 
This also naturally defines a phase space as the cotangent space P = T*Q, in 
which we can identify a pair of position and momentum coordinates {x,p) as the 
base point x E Q and p e T*Q as a hnear form on the tangent (velocity) vectors 
at X. 

On this phase space in turn is defined a Poisson bracket structure, which maps 
pairs of functions on P to new functions on P, defined in local coordinates {x, p) 
as 

A Hamiltonian, defined as a function on the phase space, then describes the 
classical evolution of the system by the flow of the Hamiltonian vector held: 

Now quantization amounts to the following. We want to map classical observables 
to quantum operators. These operators must act on wave functions in a vector 
space with inner product, that is, a Hilbert space. In the following, we will 
construct the picture bottom up; this is contrary to the way one would construct 
a quantum theory in practice, but allows us to specify the mathematics involved 
more rigorously. 

We consider a Hilbert space 5' of wave functions mapping from the configuration 
space Q to the complex numbers. The inner product on this Hilbert space defines 
normalization of wave functions. Normally this Hilbert function space is chosen 
to be the square integrable functions L'^{Q; C). 

Within this Hilbert space, wc consider a finite dimensional Lie group G within 
the set of unitary operators U : S ^ S. Such a Lie group G has a Lie algebra 
Q associated with it, which carries the standard anti-symmetric multiplication 
structure, called the Lie brackets, such that for all f,g,h e G- 

[f,9]^-[9,f] 
[f + g,h]^[f,h] + [g,h] (5.2) 

0^[fA9,h]] + [g,[h,f]] + [h,[f,9]] 

We are then looking for a set of classical observables, that is, functions on the 
phase space P, a Lie group and its corresponding algebra Q and a mapping from 
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these classical observables Oi to Lie algebra elements Oi & Q, for which the 
Poisson brackets are mapped to —i times^ the Lie brackets. In this mapping of 
the classical operator Poisson algebra to a Lie algebra, one normally also includes 
a 'quantum deformation'. This means that these relations are modified by adding 
a factor like h to introduce a quantum scale. 

This process of canonical quantization is by no means a straightforward proce- 
dure: the choice of Hilbert space, classical operators and mapping to quantum 
operators is not prescribed. One must just try and see whether the choice made 
gives a consistent scheme. 

In ordinary quantum mechanics on M", one promotes x and p to operators, which 
yields the well-known commutation relations 



Here we see that a quantum deformation of a factor h has been added. 

This is then combined with the Hilbert space S = L^(M"'; C) of squarc-intcgrablc 
wave functions over the positions Xi G M". The operators Xi and pj arc realized 
on this Hilbert space in the position representation as the self-adjoint operators 



In the outline above, we have not yet talked about constraints. These do not show 
up in the canonical quantization of a particle: there all degrees of freedom are 
truly physical degrees of freedom. In the polygon model we have to accommodate 
constraints. This can be done in essentially two ways: we can either implement 
those at the classical level and then quantize, or we can first quantize the theory 
and only then implement the constraints quantum mechanically. 

Implementing them at the classical level means that we solve the constraints clas- 
sically by constructing a fully reduced phase space. Thus we have to reparametrize 
the phase space coordinates of the system in such a way that the constraints are 
fulfilled by construction. 

We must be careful in solving the constraints: the theory was formulated as a 
generally covariant theory, which means that if we solve all constraints, we are left 

^This depends on conventions: physicists tend to add the —i, with the effect of the operators 
becoming Hermitian, thus having real eigenvalues. The operators in the Lie algebra however 
have to be anti-Hermitian to give unitary operators when exponentiated. 



[xi,pj] = ih{xi , Pj} = ihSijl. 



(5.3) 



Xi 4){x) = Xi 4){x) 




(5.4) 



5.2 Constraints in quantization 
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with a so called 'frozen time' picture. The Hamiltonian itself was a constraint in 
this formulation, so time evolution is a (semi) gauge transformation. This means 
that if we naively solve the Hamiltonian constraint, then in the fully reduced 
phase space, the Hamiltonian constraint is just a constant, so the Hamiltonian 
equations of motion become trivial. The state space of the system is then exactly 
parametrized by a complete set of constants of motion. 

This problem of a frozen time picture can be overcome by breaking general 
covariance. We can introduce a time again (often referred to as 'internal time') by 
choosing a classical observable that is monotonic in the original time parameter. 
The simplest choice for such a parameter is one of the edge lengths. Let us label 
this length and its conjugate boost by (/o,?7o)- The action 

){t)] = / dt\2r]oio + J2'^mh-uH-Y,{vfCi + WfCi] 

(5.5) 

then yields the correct equations of motion. We see that the Hamiltonian and all 
constraints have Lagrange multipliers {u,Vf G M, if/ G C) associated with them; 
the equations of motion for these multipliers result in the constraint equations, 
cf. 1(3. 9)| . If we now solve the equation H{r]o,rii) = for r/o, we obtain ?7o(^j) 
and plugging this in the action, while substituting integration over t with Iq, we 
obtain 

Sre4ih,V^){lo)] = j<^k E ^ ^ + ^ (^^) " ^ E (^/ + ^'Z C/) 

(5.6) 

We see that we have removed one pair of canonical coordinates, while solving the 
Hamiltonian constraint. The Hamiltonian is now given by i/red = — 2?7o together 
with the remaining constraints and their Lagrange multipliers. Also, some of the 
complex constraints will explicitly depend on the new time parameter Iq. 

In theory, this method of introducing an internal time will thus yield a system 
with non-trivial equations of motions. The explicit equations of motion will 
depend on the solution of f]o{rii) from H = 0, which is in practice hard to solve. 

The other way of implementing the constraints at the quantum level comes down 
to first quantizing the complete phase space and then constructing quantum 
mechanical operators C for each constraint C according to Dirac's prescription. 
A physical wave function (p must then satisfy 

C(f) = 0. (5.7) 

Again, we are faced with the Hamiltonian constraint. If this constraint is enforced 
at the quantum level, we have H (p = 0, and thus that states are constant in 



time. If we first solve the Hamiltonian constraint as in |(5.6)| , we have a highly 



complicated system, where a lot of symmetry has been removed. 
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5.3 Quantization of the polygon model 



Given the classical formulation of the polygon model in chapter ^, we want to 
find a consistent choice of quantization. 

The most straightforward choice as proposed by 't Hooft in p, is to promote the 
canonical coordinates k, rjj to operators, as in standard quantum mechanics. The 
Poisson brackets of these variables are the same as in |(5.3)| (up to a factor |). 

We now run into a problem: the Stone- Von Neumann uniqueness theorem tells us 



that the canonical commutation relations |(5.3)| have one unique representation, 
when X and p are self-adjoint operators on a separable^ Hilbert space (see e.g. 
1^, p. 65]). This representation is the Heisenberg representation and implies that 



the eigenvalue spectrum of both the /j and r]j would consist of the whole real line. 
The spectra of these operators in the quantum theory do not correctly reflect 
some properties of the original classical system, namely, the inequalities for the 
lengths k > and the triangle inequalities |(4.5)| for the boosts. We will now try 
to analyze whether these apparent problems can be remedied. 



5.3.1 Transitions 

The inequalities that lengths should be positive can be directly related to graph 
transitions. In the classical case, when a length would become negative, we 
demand that a transition takes place instead. This is because, at least classically, 
introducing negative lengths might give rise to all sorts of strange things like 
negative spatial volume or overlapping coordinate charts. These overlapping 
coordinates would also show up when we ignore transition d in figure ^.6| ; a case 
that is not covered by the inequahties li> 0. 

The question is now how we want to treat this in a quantum theory. Allow- 
ing negative lengths and ignoring transitions will clearly remove the non-trivial 
dynamics from the model. Therefore it seems that we should also restrict the 
quantum mechanical system to these length inequalities and probably rule out 
overlapping coordinates like in figure |4.4|, too. 



If we try to implement transitions in the quantum theory, we face the fact that 
transitions change the graph structure and the coordinates of the phase space. 
In the case of a polygon merging or splitting transition, they even change the 
dimension of phase space. This poses a significant problem for the formulation of 



topological space is separable when it has a dense subset that is countable. A Hilbert 
space of countable dimension is separable. If separability is not demanded, then other rep- 
resentations can be constructed, see [Q. In these representations either a; or p will not be a 
well-defined operator, although the exponentiated form will be well-defined. The spectrum of 
the other operator will still be the whole real line though. 
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a quantized theory: how should these different phase spaces be brought together 
in one quantum theory? 

In 't Hooft suggests to make an analytical continuation of the wave functions, 
such that they are continuous across transitions. Implementing the constraints on 
the quantum level should then also ensure that the wave functions are constant 
under gauge transformations. 

There are however not only the continuous symmetries generated by the con- 
straints, but there is also a discrete overcounting of physical states: one and the 
same universe can be represented by multiple completely different (on a graph 
structure level) tessellations. It might be possible that these different graphs can 
be transformed into one another by finite gauge transformations, but this is not 
a priori clear. 

This raises questions like whether one wants to keep those overcounted graphs 
on the quantum level or not and how to deal with either option. If one wants to 
remove this symmetry, then what should be the choice of gauge fixing? If one 
wants to keep these discrete symmetries, how do we construct one global Hilbert 
space and a proper measure on it? 

For the first option, a possible way to partially solve the problem of gauge fixing 
is offered in There it is conjectured that every element in the physical phase 
space can be constructed as an OPT. This would allow one to write all multi- 
polygon tessellations as OPT's by a suitable gauge transformation and thus keep 
the phase space of fixed dimensions. Then one would either need to explicitly 
find this 'suitable gauge transformation', which seems a non-trivial task, or one 
would have to show by other means that the phase space is complete. 

5.3.2 The triangle inequalities 

In the construction of the polygon model, we imposed that there is no cosmologi- 
cal constant or matter, which implies that space-time is fiat. These are sufficient 
conditions to obtain the triangle inequalities. 

As these inequalities are inherently connected to the way the model is constructed, 
it seems that they should also be imposed on the quantum level; otherwise, 
the quantum theory would have regions that are classically forbidden and it is 
therefore difficult to see how a correct classical limit could emerge. Furthermore 
the classical Hamiltonian is only defined for those boosts that satisfy the triangle 
inequalities. 

Each triangle inequality in itself does not restrict a single boost rji to a certain 
domain, but the inequalities do couple boosts at a vertex: choosing values for 
two boosts at a vertex restricts the value of the third boost to a bounded domain. 
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Even stronger, these inequalities couple between all boosts due to the fact that 
each edge is connected to two vertices and all edges to each other through 
the graph structure of the tessellation. This is in conflict with the canonical 



commutation relations (5.3), because by the Stone- Von Neumann theorem, each 



boost r]i should have spectrum M independently. This problem could be overcome, 
if we can find a symplectic transformation^ of the phase space that decouples the 
inequalities for the new coordinates. 



5.3.3 Other choices of quantum variables 

As laid out above, the coordinates k, rjj are both restricted by inequalities at the 
classical level, which gives problems when trying to promote them to quantum 
operators. A possible way to avoid this problem is to choose a different set of 
coordinates to be promoted to quantum operators. 

We can try to reparametrize our coordinates in such a way that the inequality 
constraints simplify and/or decouple; or we might be able to find a reparametriza- 
tion, which exactly parametrizes the part of phase space where the inequalities 
are satisfied. 

These methods of choosing a different coordinatization of phase space can be 
extended to a more general way to attack the problem of these inequalities. 
We can choose a set of phase space functions and promote these together with 
a quantum deformation of their Poisson algebra to operators in a Lie algebra. 
This is not a fundamentally different approach from the first; yet it does add 
the possibility of choosing an enlarged set of phase space functions that are not 
independent. 

We have tried to find a reparametrization of the boosts. The restricted pairs of 
boost triples are shown in figure |4.8| . As a way to parametrize only the triples 
satisfying the triangle inequality, we can first extract a common scale factor s 
from a boost triple. Then a triangle remains to be parametrized. We can exploit 
its symmetry by choosing three coordinate axes as in figure |5.1| . Now a point in 
this triangle can be written as x = Yl^i=i but there is a redundancy in this 
description, which can be lifted by imposing the constraint C = 1 — Yl^=i ~ 0. 
We also have to restrict the range of the coefficients: Cj G M+ suffices, because 
it implies Cj G [0, 1] together with the constraint. Furthermore, a consistent 
quantization exists on L^(]R_|_), see 0. 

This approach does not solve all problems. The new variables are a linear 
combination of the boosts, so the Poisson brackets will still be of the same form. 

•^A symplectic transformation is a transformation on a space, which has a symplectic form 
defined on it. That is in our case: a phase space with a Poisson bracket, where the symplectic 
transformation by definition preserves the symplectic form (thus the Poisson bracket). 
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Figure 5.1: Coordinates used to parametrize a triangle. 



Since the new coordinates q have a restricted domain, the same problem will arise 
from the Stone- Von Neumann theorem. Next to that, this way of parametrizing 
the boosts introduces new variables q per vertex, while each boost is connected 
to two vertices. Thus, one would have to introduce constraints of the form cj ~ 
for each boost, where 1 and 2 label the two vertices a boost is connected to. 

Some remarks can be made about an approach in the way mentioned above. We 
have seen that quantization using variables that obey the canonical commutation 
relations gives problems with the eigenvalue spectrum. There is also the prob- 
lem that the constraints (and Hamiltonian) are very complex non-polynomial 
expressions in terms of the boosts. This will probably give rise to severe operator 
ordering ambiguities when we try to implement these constraints at the quantum 
level, if suitable quantum representations can be found at all. 

To circumvent these problems, we can look for a set of phase space functions 
that do not obey the canonical commutation relations and in terms of which 
the constraints are functionally simpler. As an ansatz we propose the functions 
li, aj = sinh.{'r]j) and '~fj = cosh.{rij). These have the advantage that the linearity of 
the constraints in the lengths is preserved and the expressions obviously become 
simpler in terms of o'j,'~fj. The angles are given in terms of inverse (co)sines of 
rational functions of the <Jj,jj and thus not polynomial yet. These functions obey 
the Poisson algebra 

{ /j , (Tj } = 2 ^ij Jj 5 

{k,lj} = (5.8) 
{ > 7j } = 0, 

which are not the canonical commutation relations anymore. We now want to 
find a representation of this algebra such that / and 7 have positive spectrum. 
We have not worked out whether choosing li,aj,^j as quantum operators yields 
a consistent quantization scheme though. 
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A different approach would be to only use variables associated with a set of 
fundamental loops of the spatial manifold S^. With each such loop we can 
associate a Poincare transformation, determined by 6 parameters and there are 
2 g—1 independent loops. If we could rewrite the polygon model in terms of these 
variables, there would only be 3 constraints left. The 6(2 ^f — l) variables are 6 
more than the 12{g — l) dimensions of the physical phase space; 3 of these are 
constraints and the other 3 are gauge freedoms associated to these constraints. 



5.4 Conclusions 

The original question whether space is quantized or time is discretized seems to 
depend on the way one quantizes the theory. This fact is well known, see e.g. 
or 0, chap. 13]. 

To draw conclusions about discreteness, one must first have a consistent quan- 
tization scheme. In this chapter we saw, that the obvious choice of canonical 
commutation relations has some serious problems which cannot be overcome 
easily. Most notably the problem that the spectrum of the basic quantum 
operators conflicts with the classical length and boost constraints. Therefore 
it seems that one should not draw conclusions about the spectra of length and 
time, without first showing that a consistent quantization scheme can (at least 
theoretically) be found. 
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The constraint algebra 



In this chapter, we will do some explicit calculations to determine some properties 
of the algebra of Poisson brackets of the constraints. Next, we try to interpret 
these constraints and the transformations they generate. 

The algebra of constraints expresses how a constraint transforms under an in- 
finitesimal gauge transformation with respect to another constraint. In the case 
of first class constraints, the result should be equal to zero on the constraint 
surface (a so-called 'weak equality'), because a classical solution must satisfy 
these constraints and a gauge transformation preserves the solution and thus the 
constraint too. 

The calculation of the constraint algebra is interesting from different perspectives. 
It allows us to check that these weak equalities hold, which is a consistency 
check for the polygon model. Furthermore, if one wants to quantize the model 
and implement the constraints in the quantum theory, then this Poisson algebra 
should have an analogue in terms of operator commutators in the quantum theory. 
If they are not equal, this could give rise to quantum anomalies p. 279]. 

The constraints do not form a true Poisson algebra, because the right-hand sides 
of the Poisson brackets turn out to be non-linear expressions in terms of the 
original constraints and thus the linear span of the constraints does not form a 
closed space under the Poisson brackets. A true Poisson algebra is a vector space 
X with multiplication and a bi-linear and anti-symmetric Poisson bracket 

{.,.}: X X X ^ X, (6.1) 

satisfying the Jacobi identity. Thus when we choose a basis Cj of X, we can write 

{ei,e,} = /^.efe, (6.2) 
where the fjj are the so-called structure constants. 

Our algebra of constraints is not of this form, but can (in principle) be rewritten 
in a form where the structure constants f^j are replaced by structure functions 
fij{l,r]) on phase space. For a proof, see P, p. 8 and appendix l.A]: the Poisson 
brackets turn out to vanish on the constraint surface and can thus be written 
as linear combinations of the constraints, with coefficients given by structure 
functions. These structure functions are not explicitly given however and might 
be very complicated. 
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6.1 Useful formulas 



Before starting the calculation, we summarize the constraints and some useful 
relations, which will be used to calculate the Poisson brackets. For each polygon 
we have two constraints: the angular constraint Ce and the complex (closure) 
constraint Cz- They are |(4.7) : 



N 



ci 



^(tt - aj) - 2 71, 



Ci 



N 



Ij exp{i9j) 



where 



a,;. 



The label A = 1 . . . P labels the different polygons. The complex constraint takes 
values in C, so it is actually two real constraints. Instead of working with the 
real and imaginary part, and will be used to calculate all independent 
Poisson brackets. 



Next we have the vertex relations (4.4) , from which the following identities can 
be derived: 



dag 

dT]a 

dan 



a 



(7b 7c + Ca Oc 76 



(6.3) 



(6.4) 



dri^ Sa (Tb 0"c Sa (Tb CTc 

where the indices a, 6, c label the three different edges and angles at a vertex 



according to figure 2A_. Note that these results differ by a factor 2 from those in 
(appearing from chain-rule differentiation of the 2ri arguments). 



The rightmost expression of |(6.4)| contains only quantities that can be expressed 
in terms of boosts and angles of a single polygon at that vertex. This can be used 
to derive another useful relation between neighboring edges and angles along a 
polygon: 



fc+i 



2 (cfe+i ±isk+i) 



doik+i ^ ^ . 7fc+i 

7fc + Cfc+i C7fc 7fc+i =F i Sk+i (Tk 7fc+i 



9ak+i , ^ . 7fc 

~^ ± 2« — 

orjk CFk 



Sk+l O'k CTk+l 

Tfc 



(6.5) 
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The relation e**^'^ = —e is used throughout the calculations; together with 
the formula above, we can use it to simplify expressions a lot, as we will see. 

As most of the work needed to calculate the different Poisson brackets is very 
much alike, we will first derive two general expressions, which will be used in 
further calculations. Both expressions are of the form 



N 

E 

k=l 



(6.6) 



where we have two polygons and differentiate angles of the first (indices j) with 
respect to boosts of the second (indices k). 9j,^k denote the angles of edges as 
(4.6)| . This expression is calculated both for the two polygons being equal and 



m 



different. In these calculations we make some assumptions |(6.11)| , which can be 
found in the next section on page 

First for j, k indexing the same polygon, we have 



N 

E 



k=l 



dr]k 



N 



da. 



k=l 



-E' 



1=1 



dai 
drji 



+ e 



±i{TT-ai) 



7-1 



dai 
drji 



using the fact that angles only depend on neighboring boosts. Notice the extra 
term e^*^'^ ■''•^ : it is present because of the jump in complex angle factor, where 
the indices are taken cyclic modulo N (but not sums, as noted in |(2.5)| on page 0). 
This extra factor is 1 on the constraint surface, but non-trivial away from the 
constraint surface. Now we plug in the identity = 2i (^ — ^) and shift the 
index of one of the terms to obtain 



E' 

1=1 



dai li-i 



dr]i 



i-i 



CTl-l 



+ e 



±i(Tr-ai) 



±2 z ( e^''^ ^ - e^''- — ) + (1 



a 



aj aN ' 



+10], 



dai 



N 



dai ll 
drji ai 



doii . In 
driN aN 



(6.7) 
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In the second last step, we used relation |(6.5)| to let all terms in the summation 
cancel against each other, except for the boundary terms. 

When the j, k index different polygons, we get (again under the assumptions |(6.11)|) 
only a contribution from the common edge denoted by jc, kc- This reduces the 
number of terms in the summation: 

N 



k=l 



89, 



E 

k£{kc-l,kc,kc+l} 



-E 

1=1 



dai 



—e 



fce{-i,o,i} 



Ol^kc+k OUkc+k 



+ 



du. 



3±i(/3fc^+$M<5fec,i) 



duk,-i 



5. 



>+ie[i,i] 



2e^ 



+ 



Ckc O'kc-l 



•^ic+ielij] 



where the term e^**" = 1 on the constraint surface and we split off (1— e^**^ 
giving 



2e 



±^ '^ic6[ij] =F « — — (^j^+iefij] 



_|_ _ gTi-I-iVi^cAf^ gTiftc+i 
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±2 1 e^"^" 
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±i<I>M j 
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1 



) 



(6.8) 
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Most of the above calculation is a straightforward expansion of the definitions. 
However one must keep good track of the indices and cyclic boundary conditions, 
which make the result non-trivial: if j = N and k ^ {1,M} then the result 
reduces to zero identically. Thus we see here already that the result depends on 
the choice of labelling of the polygons from the Sk^^i and 6k^^M terms. 



6.2 Calculation of the algebra 

We want to calculate the complete Poisson algebra. All different, independent 
Poisson brackets between the constraints are given by 

{^7^cf}, {ct,cf}, {c^cf}, (6.9) 

where A, B independently run over all polygons. All other combinations can 
be obtained by using anti-symmetry of the Poisson brackets or by complex 
conjugation, for example: 



To start with, there are some constraints which are trivially equal to zero. First 
of all, each constraint commutes with itself (when also the polygon labels are 
identical: A = B). Secondly 

{C^,C^}=0 (6.10) 

because there is no dependence on the edge lengths in either constraint. 

Next we will calculate the Poisson brackets between all remaining constraints 
under some non-generic assumptions. Similar calculations were made previously 
T8[] . Here, we will perform the calculation more explicitly and establish several 



m 



corrections to the results obtained in 18 . These results we will then use as a 



starting point for calculating the general form of the algebra, thus lifting the 
assumptions |(6.11^ 



The assumptions are as follows: 

no single polygon has two edges that have to be identified and (6.11a) 
no two different polygons have more than one edge in common. (6.11b) 

These assumptions are a theoretical, but also practical restriction, as they do not 



cover the 'one polygon tessellations' (OPT's, as in |13|). There every edge of the 
polygon has to be identified with another edge of the same, single polygon. OPT's 
are the simplest tessellations possible and therefore a practical set of tessellations 
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Figure 6.1: Label scheme for the calculation of the constraint algebra 



to parametrize phase space with. If we want to explicitly use these OPT's as a 
starting point for a quantum theory, we therefore need the full Poisson structure 
without the assumptions. 

We will calculate the explicit expressions for the general Poisson structure without 
these assumptions and show that the constraints do also close in that case. This 
will be treated in section 0| and a summary of all results will be presented in 



section EA 



When we calculate brackets between two different polygons (labeled I and II), 
the result will be zero unless these polygons have an edge in common. This is 
because the constraints only depend on variables locally at the border of the 
polygon. Hence with our assumptions we only have to look at cases where two 
different polygons have exactly one edge in common. 

In the calculations that follow, we will parametrize the two polygons as in fig- 
ure BTTI. We label lengths, boosts and angles of polygon I with /, rj, a respectively 



and j indexes these. For polygon II we use m, (3 with indices k. The indices 
of the edges of polygon I and II that join each other will be denoted by jc 
and kc respectively. Labelling will always be counterclockwise oriented, see also 
chapter ||. 

We start by calculating the Poisson bracket of the angular constraint with the 
complex constraint of the same polygon. As there is no dependence on the lengths 
in the Cg, we only have to differentiate it with respect to the boosts and thus Cz 
only with respect to the lengths. This yields 

N 



N 



k=l 
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0. 



(6.12) 



In the last step we used equation |(6.5)| to rewrite the result to the same form 
as in [0. It coincides with the latter expression up to a factor 2e*^'^~°^^; the 
exponential is due to a different definition of the complex constraint and the 
factor 2 probably due to a calculation error, as noted below |(6.4)| . 

When we take the complex constraint together with its complex conjugate, we 
get terms from both differentiating the first with respect to the boosts and the 
second with respect to the lengths and vice versa: 
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using 1(6.7) 



(6.13) 



This result is somewhat different than that in [0]: in the first term we have a 
relative minus sign instead of plus and the second "boundary" term is added; it 
arises due to the cyclicity of the edge and angle indices, but one has to be careful 



with the interpretation of sums of angles as noted in |(2.5' 
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Notice that the correctness of the minus sign is easily checked by complex con- 
jugating the complete Poisson brackets, because is should satisfy 



{ c^. , a } = { a , } = - { , a } . (6.14) 

We can check that the second term does also obey this condition. 

Next we have the Poisson brackets between two different polygons. Using equa- 
tion [6.8)1 we get 

M 

9n ,^mk e 

de 



k=l 

M 



-IE 



k=l 

1 



k 



-I e 



= (1 - e**") + e-*^i 5k^,M) ~ (6.15) 

by straightforward calculation, where we have to remember that the indices are 
cyclic and thus trivially ^jcG[i,Af] ~ ^jc+i&[i,N] = 1 — 1 = 0. Comparing again with 
0, we have found an additional boundary term arising from cyclicity of the 
indices along the polygon. 
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The calculation of Poisson brackets between two complex constraints gives 

N M 

J2 h > E 
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3,k 
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k=l 



) 
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inserting |(6.8) 
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where we have defined 
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0, 



(6.16) 



AT 



Zt = Y.^3^'' (6.17) 

j=fc 

as the partial complex constraint of polygon A from edge k. Also, we used the 
relations 

'^ice[ij] - (^jc+iefij] = "^jjc - ^iciv and '5j,e[i,i] = ^m<^.^\- (6-18) 
Completely analogously, we find for a complex and complex-conjugate constraint 

= H Sj^^N 6k,,M 



3c 



+ 



+ 



[zl^,5,^,, + e-^'^^ZlS,^,,, 



M 



(6.19) 
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Figure 6.2: Rotation of the plane of a polygon and its action on edge lengths and 
extrinsic curvature surrounding the polygon. 



Both of these Poisson brackets yield expressions that vanish on the constraint 
surface, but do not identically vanish on the whole phase space, unlike the results 
in [|13]: when one of the indices jc, kc is equal to the first or last of its polygon. 



the result is non-zero. 

As already pointed out in [0, the Hamiltonian flow generated by the complex con- 
straints should (on shell) be related to Lorentz transformations of the coordinate 
frames of the polygons. This is a residual gauge freedom that is left after fixing 



the constraints in |(3.8)| . With this interpretation, one would not expect the 
Poisson brackets between these constraints to vanish identically (as is the result 
1T8|): the complex constraint should be related to the constraints and their 



m 



corresponding gauge freedom to the Lorentz transformations of the local frame 
of the polygon. 

Suppose now that the complex constraint in polygon I does not close by an 
amount d in the ?/-direction (C^ = id ^ 0) and we generate a Lorentz transfor- 
mation in polygon I. This Lorentz transformation can be seen as a "space-time 
rotation" of the plane of polygon I as in figure |6.2| : the (infinitesimal) boost of 
the Lorentz transformation corresponds to an angle that the plane of polygon I is 
rotated by. This will affect the lengths and boosts of all edges along that polygon 
and also the lengths of edges connected to those. We now consider the effect of 
this transformation on a neighboring polygon II. We suppose that they have an 
edge in common, which has indices jc, kc in polygons I and II respectively, where 
we choose jc = N,kc = M. The implicit assumption is made that the closure 
deficit Cl is located between edges j = N and j = 1. Now the position of vertex 
1 (which is located between edges j = N and j = 1) depends on whether we 
view it as the first or last point of the constraint C], see figure |6]^, where we 
have implicit complex coordinates z = x + iy and the origin is placed at vertex 
1 viewed as starting point of the complex constraint. 

This will give rise to a different change of the intersection line under this Lorentz 
transformation and thus a change in Cl^. Figure |6]^ shows the geometrical 
picture in the {y, t)-plane. The lines are the planes of polygons I and II and 
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II t' 




Figure 6.3: Action of a Lorentz trans- 
formation on a complex constraint, (y, t) 
and {y', t') denote the coordinates within 
polygons I and II respectively. 




Figure 6.4: Polygon with non-closing 
complex constraint: vertex 1 shows up 
on both sides of the deficit 7^ 0. 



their intersection point is the common edge. Now vertex 1 shows up as this 
intersection point, but also at a distance d apart. 

The action of an infinitesimal Lorentz transformation of polygon I by an amount 
diy is shown as the dashed lines. The upper one is rotated around the origin of 
polygon I (which is vertex 1) and the lower one around vertex 1 as viewed after 
going around polygon I, after which it is shifted by a distance d. These generate 
different new intersection points between polygons I and II. We can now calculate 
the distance between those intersection points explicitly, which yields 

smh(/ 7]) 

This matches the result 

{C'I,Cl} = '^Cl5,^,^='^, (6.21) 

where we have inserted that all constraints close, except C\ = id and thus all 
other terms from |(6.16)| vanish and = 2 tt. 



Thus this geometrical viewpoint seems to indicate that indeed the Poisson bracket 
{ C\ , C^^ } is not identically zero outside the constraint surface, if we interpret 
the complex constraints as generators of Lorentz transformations. 
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6.3 Generalizations 



We can generalize the assumptions |(6.11)| made in the previous section to the 
arbitrary case, which will recover the explicit full Poisson structure. We will 
take two different approaches: the first one is to consider the Poisson brackets 
of the constraints with the Hamiltonian, using the Hamiltonian flow in |(4.8 



The Hamiltonian is the sum of the angular constraints, so for an OPT we have 
H = Cg. This will then give us the Poisson bracket {Cg , C^} iia the case of an 
OPT without any assumptions made. The second approach will be to lift the 
assumptions in a general setting. The results of these approaches are then shown 
to match in the case of an OPT. 



6.3.1 The case of an OPT 



To calculate the Poisson bracket {H , Cz}, we first need a formula relating the 
growth of edges along a polygon. This formula and its derivation is analogous 
to 1(6. 8)| . However we note that in the definition of the a, 7 in |(6.22)| , the factor 



2 is missing with respect to our original definition |(2.4)| , so cTj = sinh(r7j),7j 
cosh(77i). We have 
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(6.22) 



We can use the previous relation to rewrite {Cz , H} (the time evolution of the 
complex constraint) in the same way as we did with |(6.12)| : 

{Cz,H} 




N 



N 
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= e*^i - 2 (1 - e'^^) using (622) 

^r.l^fLtil + ^^^Ve^-lj^O (6.23) 

Here we recognize the left part of the final expression as minus the velocity of 
vertex 1 in the coordinate-frame of the polygon (in complex notation, as was used 
for the formulation of the complex constraint). The factor e*^^ — 1 precisely gives 
the difference of viewing this point as the first or the last when going around the 
polygon. Thus this intuitively matches the expected result of the time evolution 
of the closure constraint. 



6.3.2 The case of a general tessellation 

We now want to generalize the calculated Poisson brackets between two polygons 
to the arbitrary case without assumptions |(6.11) 



The underlying idea is as follows: we can split the Poisson bracket expression into 
parts (see |(6.25]| for a precise specification), taking account for each common edge 
separately. For this, we first have to show that there are some restrictions to the 
ways two polygons can have multiple edges in common, such that we can make 
such a splitting. For this to work, different common edges should not interfere 
with each other in the calculation of the Poisson brackets. 

To show that the "condition" above is fullfilled, let us consider two polygons as 



in figure 3^, which have at least an edge in common for edge indices jc, kc- The 
polygons X, Y and Z can be copies of polygons I or II and as such can give rise 
to multiple common edges. 

The only non-vanishing terms in the Poisson brackets appear from the differen- 
tiation of one term in a complex (conjugate) constraint (say polygon II) to the 
edge length and then the other constraint (polygon I) differentiated with respect 
to the corresponding boost. This means differentiation of an angle of polygon I, 
as the angles are the only quantities in the constraints that depend on boosts. 

If two common edges cannot appear in such a way that their 'split' Poisson 
brackets would both contain a specific angle term of polygon I being differentiated 
to a boost of the same edge of polygon II, we can replace the Poisson brackets by 
two separate Poisson brackets where we treat the common edges one at a time. 
Note that here an 'edge' is a boundary element of a polygon, so each pair (kjiji) 
is assigned to two edges, when the edges are viewed as boundaries of polygons. 
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Thus the "condition" above can be stated in a more precise way as follows: 

each edge from polygon II has at most one way in which it is connected 

to each angle of polygon I (and the converse, which follows from (6.24) 

interchanging polygons I and II). 

The proof of this statement hinges on the fact that the graph of the tessellation 
has trivalent vertices. We will first give an insightful proof followed by a more 
detailed proof. 

The statement can be seen to hold in the following way. If an angle would 
be connected twice to an edge, then this edge has its endpoints on one vertex. 
Because of trivalency of the vertices, both ends of that edge must be neighboring 
at that vertex and thus the edge encloses a polygon in itself. This is clearly a 
contradiction. 

For the detailed proof, we start off with the case of an angle of polygon I and an 
edge of polygon II. Proving |(6.24)| boils down to a tedious but straightforward 
treatment of all possible scenario's. With some symmetry arguments, we can 
however reduce this a bit. Suppose that an angle of polygon I is connected twice 
to an edge of polygon II. It must be an angle at one side of a common edge. 
Because of symmetry we can take angle a as shown in figure |6.5| without loss of 
generality. 

Now edge kc cannot be connected to a more than once, because it is already 
joined with edge jc and thus cannot be joined with — 1 and it cannot be equal 
to edge kc + I. Also edge kc + I cannot be connected more than once, because 
if it would be connected to jc — 1 (as b), then kc would be connected to fee + 1 
(as a). Finally, edges a and b obviously cannot be equal, so condition |(6.24)| is 
fuUfilled for angles of polygon I. 

To apply the idea as outlined above to the generalization of the Poisson brackets, 
we first have to explicitly define what we mean by a splitting of a Poisson bracket 



kr-1 




II 



Jc + l 




Figure 6.5: Multiple common edges between two polygons. 
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into separate parts for each common edge. Let us denote the pairs of common 
edges of polygon A and B by {jc, kc)c=i..nc- We now introduce the notation 

{C^.C^], and {C^,C^}^. (6.25) 

The left expression shall denote the Poisson brackets of two constraints, where no 
edges are identified. This expression will thus only be non-zero when A = B and 
then yields the 'internal terms'. The right expression shall denote the Poisson 
brackets of these constraints, where no internal contributions are counted (so we 
treat A B) and only the edge pair [jc, kc) is treated as common; the others 
are treated as not being common. Thus by their definitions, these expressions 
are equal to the corresponding expression calculated in section |6]^ under the 
assumptions |(6.11)| : { , for the case A = B and { , }^ for the case 



A B. We will generalize the Poisson brackets already calculated to arbitrary 
conditions, by proving that one can split the general Poisson brackets as 

rtc 

{C\C^} = {C^C^}^ + Y,{C\C^}^. (6.26) 

c=l 

Given two polygons I and II, we denote their pairs of common edges by {jc, kc)c=i..nc 
and first investigate the Poisson brackets { 9^^ , C^^ } : 

rtc 
c=l 

c=l 

i'=i 

= {e], Cf } (6.27) 

Working backward from the desired result, we first only write down those terms 
that will contribute for each common edge. In the next step, we use the fact that 
each combination of angle a of polygon I and edge m of polygon II does show 
up not more than once in the full Poisson bracket, as stated in |(6.24)| . It also 



does show up exactly once, as it is part of a single common edge bracket. All 
remaining terms in | , Cl^ } are zero now, because we have exactly taken into 
account all terms contributing to all common edge Poisson brackets and there 
are no internal contributions, because we considered two different polygons. 
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In the case that both constraints are from the same polygon, we can essentially 
repeat the argument, except that we have to add the terms that arise from the 
internal contributions: 



c=l 

N 

{ -Oik helij] - ttfc+i 4+ie[i,j] 1 k e*^'' } + 

k=l 

E [ { '^^■^^[i.il ' ^'^-i + ^^''^ } + 



c=l 



E { . c!' } 



} 



(6.28) 



The extra terms do not change the argument, because also with these added, 
condition (6.24) still holds and all terms that we started with are exactly those 
that will be present in the full Poisson bracket. 

Now that we have shown in |(6.27)| and |(6.28)| that a split can be made for the 
brackets | , } for all A, B, we can easily extend the result to all possible 
Poisson brackets. The apphcation of |(6.27)| to {Cj , Cl^ } is trivial by plugging 
in j = and yields 

{Clef} 

ric 



c=l 



c=l 



using 1(6.27) 



(6.29) 



We see that the sum over all joined pairs of edges (jc, ^c) reduces to only contri- 
butions from the pairs where kc = 1 or kc = N. 
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In the same way, we can generalize { , j using |(6.28)| , which gives 



c=l 



_ / (^N 7i + ci ^1 In _ .1n\ _ ^ie^x 

+ ^ (1 _ e^*«) (4^,1 + e-'P^ ~ 0. (6.30) 



c=l 



To generahze the Possion brackets of two complex (conjugate) constraints to the 
case of multiple common edges, we first rewrite one complex constraint to a sum 
of 6j by the chain rule for differentiation and obtain 

N M 

i=i fc=i 

N M 

j=l fc=l 

j=\ c=l fc=l c=l 

= E { ' }c (6-31) 

c=l 

The Poisson brackets where one of the constraints is a complex conjugate and 
when both constraints are from the same polygon, are all derived in an analogous 
way. For the explicit expression one must insert the already calculated brackets 
with assumptions. See |(6.32b)| and |(6.32c)| for the complete expressions. 
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6.3.3 Comparing for consistency 



Now we can do a consistency check with the results we have found. The Hamil- 
tonian can be written as the sum of all angular constraints. The Poisson brackets 
of this sum of all angular constraints with one complex constraint should thus 
match expression (6.23)| . Indeed we find that they do match: 



+ V (1 - e^*«) 4.1 + ^ 5,,m\ 



where we sum over all angular constraints, so all separate summations over 
common edges reduce to one summation over all edges of polygon B. This sum 
then reduces to the two terms k = 1 and k = M, which leaves us with 



(1 - e'*«) 
(1 - e'*«) 



o-Af 7i + ci (Ti 7m , .7Af , 1 e" 
- + I \ h 



Si CTl CTm 

1_ ^71 - 1 

Si 



aM siai siOm 



lM-^\ . .7a/ - 1 

+ Ci 1 + I- 



Vl + Ci Vn 

Si 



using 



7-1 



a 



inserting (6.23 



Note that we replaced ^'s by $'s when inserting |(6.30)| and vice versa when 
rewriting the result back to |(6.23)| , but these are the same as both constraints 
describe the same polygon. 
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6.4 Summary 



The results of this chapter are summarized below. The same conventions have 
been made and the sums are again over all pairs of common edges between 
polygons A and B. In the case that A and B are equal, each pair will thus 
show up twice: if the pair {jc, kc) = {x, y) exists, then also {jc, kc) = {y, x). 



A,B 



+ E ( (1 - e^*") (^^c,i + e-'^^ 4e,M) ) (6.32a) 



c=l 



c=l 



Okc 



+ 



S/Si o, 



(6.32b) 



(1 - e-'^^) + (1 - e"^^) 



c=l 



^jc,N - —— '^z ^ko,M 



Okc 



'3c 



+ 



(1- 



fee "jc,iV 



(6.32c) 
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6.5 Interpretation 



We would like to interpret the Poisson bracket relations we have found. The 
angular and complex constraints are first class constraints and should correspond 
to some residual constraints and gauge freedom left after the (partial) fixing of 
the constraints 7i, W in |(3.7)| and |(3.8) 



Intuitively, the residual constraint left from the Hamiltonian density Ti is the sum 
of deficit angles at the vertices. Superficially, this seems to match the angular 
constraints, but in these the sum of deficit angle at a vertex is distributed among 
the angular constraints of polygons that join at that vertex. 

The complex constraints should be related to the residual gauge freedom of 
Lorentz transformations of each polygon frame. With the help of the dual 
graph formulation of appendix 0, we can try to see whether this intuitive idea 
matches the explicit complex constraint definition. In the dual graph, a Lorentz 
transformation of a polygon frame / G -F amounts to shifting the corresponding 
vertex v & V m the dual graph by the action of the Lorentz transformation L 
when we choose v to be at the origin (1,0,0). 

The boosts r]i of edges along the polygon are encoded by half the arc lengths 
of the geodesies emanating from v. Thus if we make an infinitesimal Lorentz 
transformation, the change in boosts will be half the "length" of the Lorentz 
transformation projected on the direction of the outgoing edge (as seen in the 
dual graph). This means that a Lorentz transformation of boost dp in the x 
direction in the polygon frame should affect the boosts rji as follows: 

drji = cos(0) dz/, (6.33) 

with (f) the angle that edge i makes with the x axis in the dual graph. This angle 
is = 6'j + |, because the direction of the boost is perpendicular to the direction 
of the edge in the normal graph. 

We can compare this to the gauge transformation generated by the complex 
constraint. The action of that gauge transformation is given by the Hamiltonian 
flow; the Poisson brackets of a quantity with the constraint gives the action that 
the corresponding gauge transformation has on that quantity. For the boosts we 
find 



dr], = {r],,a} = -l-^ = = -i (cos(^i) + I sin(^,,)) . (6.34) 



Comparing this with |(6.33)| , we see that the imaginary part of generates 
Lorentz transformations in the x direction and the real part in the negative y 
direction, at least for the boost variables. 

Unfortunately, the action on the length variables is much more complicated. To 
compare the action on the lengths of the complex constraint with the action of a 
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Lorentz transformation, we can consider the action of a Lorentz transformation of 
the polygon frame geometrically and check whether that matches { /j , }. In the 
geometrical picture of the tessellation graph, one expects that an (infinitesimal) 
Lorentz transformation acts on the variables as a "rotation in space-time" of 
the polygon plane, see figure |6.2| . This changes the intersection lines with its 
neighboring polygons, from which the action of the Lorentz transformation on 
the lengths can be deduced. The explicit formulation of this action is also 
quite complicated and we have not investigated whether this geometrical picture 
matches the Poisson brackets. 
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Conclusions 



We have looked at the polygon model as a way to describe 2 + 1-dimensional 
gravity. This model makes use of the local flatness of 2 + ID space-time in 
a remarkably efficient way: the polygon model can be formulated as a finite 
dimensional system with a completely explicit description of its evolution in terms 
of growth of edge lengths and polygons undergoing transitions at a discrete set of 
times. This makes the model very useful for exploring features of gravity already 
in the classical context, e.g. by simulation. 

In well known other descriptions of 2+1-dimensional gravity, the system is defined 
more implicitly. For example in the ADM formulation with York time, one must 
solve the conformal factor of the metric from a elliptic partial differential equation. 
This has no known explicit solution for genus g > 1 fT^ • Witten's approach seems 
to use the frozen time picture, where the dynamics cannot be recovered from the 



system||T5|, 0. These problems are not present in the polygon model formulation. 

There are however also some drawbacks of the formulation of the polygon model. 
The polygons that together form a spatial slice, are of finite size and internally 
fiat by construction. In this construction we have solved the field equations of the 
gauge constraints. The finite dimensional system which is left, is (by definition) 
non-local and still has constraints. This yields difficulties. 

The problems of this non-locality are most clear in the Poisson bracket struc- 
ture of the constraints. The constraints for each polygon are already non-local 
expressions and the Poisson brackets of constraints between different polygons 
are complicated expressions, which do not vanish off shell. Even stronger, they 
do not form a proper Poisson algebra, because the brackets are not linear in the 
constraints. 

We have shown though, that the complete Poisson bracket structure can be cal- 
culated explicitly and does vanish on shell. This was expected, as the constraints 
should be a residue of the field constraints in the unreduced canonical formulation, 
which are first class. 

The features of the polygon model that allowed an explicit formulation in the 
classical theory become problematic when we turn to the question of quantization. 
First of all, a full description of a state depends partially on a discrete graph 
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structure, which changes under transitions. This cannot easily be encoded in a 
quantum formulation. Furthermore, there are inequality constraints which are 
coupled in a very non-local way. This is due to the fact that the Hamiltonian is 
a highly non-trivial function of the basic variables, which is only defined when 
those inequalities are satisfied. Because of these inequalities, we cannot lift the 
simple Poisson bracket structure between those coupled variables to a quantum 
theory, as that would violate these triangle inequalities. 

We have shown that several ways of attacking the problem of finding a consistent 
quantization run into one or more of the difficulties above. Thus the question 
of whether the polygon model obtains a discrete time and/or space spectrum 
after quantization seems to be premature. It not only depends on whether 
one first quantizes and then implements the constraints or vice versa, but more 
fundamentally on whether a satisfactory quantum theory can be formulated in 
either approach. If we first quantize, we face a complicated Poisson structure 
which has to be implemented in the quantum theory. If we first solve the 
constraints, we face the problem of finding an explicit global solution to those. 

Both approaches, although meeting some major obstacles, do not seem com- 
pletely hopeless. The dual graph formulation (as in appendix 0) gives a nice 
abstract interpretation of the polygon model, which might be useful to reduce 
the constraints in the classical theory in terms of holonomy-loop variables. This 
dual graph formulation depends on the absence of particles however. On the 
quantum theory side, one might consider other variables than the canonical ones. 
This leaves the question of whether this polygon model can be quantized to a 
difficult but open one. 
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The dual graph 



In this appendix we will present the dual graph of a tessellation in the polygon 
model. For this we will first introduce the concept of a dual graph in its general 
setting and also introduce hyperbolic space. It turns out that we can assign 
a meaning to the dual graph of a tessellation, such that it becomes a graph 
embedded in a hyperbolic space: a surface with constant negative curvature. 

A.l The dual of a planar graph 

The concept of a dual graph is a very general one within graph theory and has 
applications in a wide variety of fields. It does not add new information to the 
graph, but allows one to see some features more quickly by using a different 
representation of the same graph data. 

To introduce the dual graph, we start with a graph G = (y,E), where V is the 
set of its vertices and E the set of edges. Each edge is normally represented as 
a pair of vertices: the ones that it connects. If we now embed our graph in the 
plane, the dual of that graph is defined by associating a vertex to each region of 
the plane and connect two vertices with an edge when their regions in the original 
graph had an edge as a common boundary. One must however notice two things: 
firstly, not every graph can be embedded in the plane without intersecting edges, 
so this definition does not yield a dual graph for each graph (V, E). Secondly, a 
graph can have different possible embeddings in the plane, resulting in different 
dual graphs; see figures |A.1| and [A.2| for an example. 

This non-uniqueness of the dual graph can be lifted if we add faces to the graph 
information. We now write G = {V, E, F) with F a set of faces and each edge also 
containing information about the pair of faces it connects. We want this graph 
with added faces to uniquely describe a graph on some genus g surface. This 
means that we have an orientation^ in the graph from the underlying manifold 

"'^The orientation can only be defined locally in the case of a generic surface. However we 
only consider orientable surfaces, which allows us to extend this local orientation to a global 
one. 
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Figure A.l: A graph embedded in the 
plane with its dual graph (dotted lines) . 



Figure A. 2: The same graph em- 
bedded differently with a different dual 
(dotted lines). 



and that at each vertex and face, we have an ordered hst of the edges that 
"connect" to that vertex or face. 

With this added structure, we can assign a unique dual graph G to each graph 
G = (y, E, F) by mapping vertices of G to faces of G and faces of G to vertices 
of G. The edges are mapped to themselves with the pairs of vertices and faces 
they connect interchanged. So we have G = (V, E, F) = (F, E, V). 

Given this definition we can observe some facts. These only depend on the graph 
structure, because at this point we have not yet added any additional information 
to the graph, like edge lengths. 

When we take the dual of the dual graph, we get back the original graph: 

G^G. (A.l) 

From the geometrical picture of viewing the graph embedded in a surface, we see 
that a graph G and its dual G represent graphs on surfaces of the same homotopy 
class. When we look at the Euler characteristic x of the graph, we find the same 
result, because we have 

X = 2-2^ = #y-#E + #F, (A.2) 
where g is the genus of the (orientable) surface. 
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A. 2 Hyperbolic space 



Hyperbolic space is a Riemannian manifold. Its most common definition (and 
most practical here) is as follows. For n > 2, start with an (n + l)-dimensional 
Minkowski space, thus R""*"^ with the metric^ r]'^'^ = diag(— 1, 1, . . .). Now n-di- 
mensional hyperbolic space is defined as 

IT = {xe M"+^ I r]{x, x) = -l,xo> 0} (A.3) 

with the metric induced from the underlying Minkowski space. 

The space thus constructed is a true Riemannian manifold: the original metric 
was a pseudo-Riemannian metric, but its restriction to H" is positive definite. 
Furthermore H" is a space with constant negative (sectional) curvature —1: at 
any point p and any 2-dimensional plane in the tangent space of p, the scalar 
curvature is —1. 

We will be concerned specifically with the case n = 2: the hyperbolic plane. 
The hyperbolic plane has some other well-known representations, besides the 
definition [(X3l given above. One is that of the complex halfplane 

(J^X'^ -\- dip' 

{x + iy eC\y>Q], ds^ = (A.4) 



and the other is the unit disc 



{ix,y)eR'\ Ux,y)\\<l}, ds' = ^''^J \ (A.5) 



2 dx^ + dy"^ 



These two models are known as the Poincare plane and disc respectively. 

Hyperbolic space is a maximally symmetric Riemannian manifold. That means 
that its isometry group is as large as can possibly be for a Riemannian manifold, 
namely, of dimension |n(n + 1). We will further investigate this symmetry from 
the viewpoint of our original definition |(A.3)| . The symmetry group is the identity 



component SO^{n,l) of the full Lorentz group 0(n, 1). The group SO^{n,l) 
consists of the time- and space-orientation preserving Lorentz transformations on 
the Minkowski space M""*"^ and thus preserves the metric rj and the sign of Xq. 

This group SO^{n, 1) can now be seen to have a group action on H": an element 
L G SO~^{n, 1) defines a map L : R"+^ i— > that can be restricted to a map 
L : H" H", because SO'^{n, 1) precisely preserves the quantities that define 



^This can be interpreted as the metric on Minkowski space itself and the differential metric 
on the tangent space as Minkowski space is a vector space and thus isomorphic to its tangent 
space. 
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H". This map L then also naturally defines the map DL : TH" i-^ TH" on the 
tangent bundle. 

The elements L G SO^{n, 1) exactly correspond to all possible isometries of H". 
This is because the group SO~^{n, 1) has precisely |n(n + 1) degrees of freedom: 
n degrees correspond to translations of the hyperbolic plane and — 1) cor- 
respond to rotations around a fixed origin. 

Looking in detail at how SO~^{n, 1) acts on H", we see for n = 2 that pure Lorentz 
boosts L^ ylr]) perform translations in H" by a distance of rj in the y-direction 
(locally around the "origin" (1, 0, 0); because of curvature the action is somewhat 
more complicated away from the origin). The rotations Lr{a) perform a simple 
rotation around the origin. 

A. 3 Embedding the dual graph in 

We have first constructed the dual graph of a tessellation. Let us now assign 
geometrical data to this dual graph. 

We will do so by looking at the path that a Lorentz transformation maps out in 
hyperbolic space when it acts on (1, 0, 0). Each Lorentz transformation L can be 
decomposed into pure rotations and boosts in the following way: 

L = mi)LMLM2)- (A.6) 

This can be interpreted (in reverse order) as rotating the local frame over an angle 
01 to the direction of the boost, then boost with parameter rj and then rotate 
the boosted frame to the desired orientation over an angle 02- The continuous 
transformation of (1,0,0) under the one-parameter subgroup 

L{t) = L,(0i) L,(r7t) L,(02), t G [0, 1] (A.7) 

will then be a geodesic in H^; figure [A.3| shows an example of such geodesies. 

Now each edge in the dual graph corresponds to a Lorentz^ transformation L 
of one frame to the other it connects. Actually, we should consider oriented 
edges, because when we traverse the edge in opposite direction, the Lorentz 
transformation is L~^. Thus each edge can be viewed as an isometry of H^, 
which in turn can be viewed as geodesies in H^. 

A triangular polygon in the dual graph corresponds to a vertex in the original 
graph. Thus traversing its three enclosing edges means going around that vertex 
and the composition of corresponding Lorentz transformations should be the 

•^We can more generally assign to each edge a Poincare transformation, by also taking into 
account the way coordinates are shifted after crossing an edge. 
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Figure A. 3: The hyperbolic surface Figure A. 4: Angles of a hyper- 
El^ with red lines showing a triangle: bolic triangle, 
the successive action of boosts and ro- 
tations, forming a closed holonomy. 



identity. This means that, again viewing these edges as geodesies in H^, they 
form a closed triangle and we can thus naturally associate a graph structure with 
it, see figure |A.3. 



We can now try to interpret the geometrical data of this graph embedded in 
H^: the dual graph lives in a Riemannian manifold, so concepts like lengths and 
angles are well-defined. Some subtleties will show up however, due to the fact 
that these geodesies do not contain all information of the Lorentz transformations 
of the edges. Each Lorentz transformation L is fully characterized by the way 
it acts as an isometry on H^, but if we only look at its action on a single point, 
information is lost: the way it acts on other points. This can be retrieved by not 
only giving the action of L on a single point, but also specifying its action on the 
tangent space at that point. 

When we now draw a triangle in the dual graph like in figure |A.4|, we find that 



each edge has geodesic length \2r]\. The outer angles of the triangle are in direct 
correspondence to the angles ctj at a vertex in the original graph: they describe 
the rotation that we have to make, when we "walk" along a geodesic and turn to 
the direction of the next geodesic. The inner angles ctj of the triangle are related 
to these as 

ai = Tc — ai. (A. 8) 

With this identification, we see that at each vertex in the dual graph, we have a 
2 71 neighborhood as expected: 

N N 

J2ai = Y,^-^i = Ce = 2n (A.9) 

i=l i=l 
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by identification of tlie vertex in tlie dual grapli witli a polygon in the original 
graph and using the angular constraint. 

There is also an interpretation for the Hamiltonian (4.9) in the dual graph by 
considering the area of all faces in the dual graph. We start from the fact that 
the area of a triangle in hyperbolic space is determined by its inner angles oii 
only, as the deficit angle from vr: 

3 

A = Ti-^ai. (A.IO) 

i=l 

When we now sum the area of all triangles in the dual graph, we obtain 



1=1 / nev 



The boost parameters 2ri cannot be interpreted as simply the geodesic lengths: 
T] can be negative and in the geometric picture of figure [A.4| this corresponds to 
walking the geodesic line in the opposite direction. 

If we have a vertex where the signs of the three boosts are not the same (called a 
'mixed vertex' in ||13[); this implies that the interpretation of cij = tt — ctj as the 
inner angles of the triangle becomes problematic. 

A way to solve this is to consider again the Lorentz transformations corresponding 



to each edge. These can be decomposed as in |(A.6)| . In decomposing these we 
can always choose rj positive: if it is negative, we can replace the parameters 
by (01, ri, 02) ^ (01 ± TT, — ?7, 02 ± vr). In this way we get a graph with all edge 
lengths positive and the oriented sum of angles at every point a multiple of 2 tt. 
If we translate this interpretation back to the original graph, it means that we 
have reversed the orientation of the edges. This interpretation thus moves the 
signs from the boosts onto the lengths, which looks similar to the variables Xi^pj 
introduced in |p. This will affect the definition of the Hamiltonian and the 
equations of motion, but we have not analyzed further the implications. 

See also the discussion of this problem of the interpretation of mixed vertices in 
P, p. 47]. 
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